Intervalley-Scattering Induced Electron-Phonon Energy Relaxation in Many- Valley 

Semiconductors at Low Temperatures 
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We report on the effect of elastic intervalley scattering on the energy transport between electrons 
and phonons in many- valley semiconductors. We derive a general expression for the electron-phonon 
energy flow rate at the limit where elastic intervalley scattering dominates over diffusion. Electron 
heating experiments on heavily doped n-type Si samples with electron concentration in the range 
3.5 — 16.0 x 10 25 m -3 are performed at sub-1 K temperatures. We find a good agreement between 
the theory and the experiment. 
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Since the low temperature hot electron experiments by 
Roukes et al. |lj , the energy transport between electrons 
and phonons has continued to be a topical subject. Re- 
cently, there has been significant experimental and theo- 
retical interest in the electron-phonon (e-ph) energy re- 
laxation in metals and semiconductors at low tempera- 
tures 0, S S S H 0, @. The understanding of ther- 
mal e-ph coupling is important for several low temper- 
ature devices such as microbolometers, calorimeters and 
on chip refrigerators 0, . This coupling plays also an 
important role in correct interpretation of low tempera- 
ture experiments [j| and the e-ph energy relaxation rate 
gives direct information about phonon mediated electron 
dephasing 10]. 

Interaction between electrons and phonons is strongly 
affected by the disorder of the electron system and, there- 
fore, the problem is commonly divided into two special 
cases: pure and impure (or diffusive) limit of e-ph in- 
teraction. The cross-over between these two regions is 
defined as ql = 1 , where q is the phonon wavevector and 
I the electron mean free path. If the whole phonon system 
is to be considered then the phonon wavevector can be 
conveniently replaced by the thermal phonon wave vector 
qr = ksT/hv, where T is the temperature of the lattice 
and v the sound velocity. Recent theory for single- valley 
semiconductors !*| predicts that the e-ph energy relax- 
ation is strongly enhanced when the system enters from 
the pure limit (ql > 1) to the diffusive limit (ql < 1). 
The behavior is the opposite in comparison to metals 
where it is well known, since the pioneering work by A. 
B. Pippard [llj . that the disorder of the electron system 
tends to suppress the e-ph energy relaxation (see also 
Ref. |2]). In semiconductors, due to small electron den- 
sity, the e-ph interaction can be described by deformation 
potential coupling constants, which do not depend on the 
electronic variables, while in metals the coupling strongly 
depends on the electron momentum 0] • This fundamen- 
tal difference eventually leads to disorder enhancement of 



the relaxation in the diffusive limit in single-valley semi- 
conductors [8(. 

In many-valley semiconductors the situation is further 
altered due to intervalley scattering, which is the topic of 
our work. Due to lack of screening the e-ph energy flow 
rate is strongly enhanced in many valley semiconductors 
in comparison to single valley ones at diffusive low tem- 
perature limit. We approach the e-ph energy transport 
problem by first considering the phonon energy attenu- 
ation rate due to electrons (or phonon-electron energy 
relaxation rate). This procedure is attractive, because it 
enables straigth forward comparison between our work 
and previous literature, which has concentrated mainly 
on ultrasonic attenuation [H Q E E E3 • We de- 
rive expression for the total e-ph energy flow rate (by 
using the phonon energy attenuation rate) and perform 
low temperature electron heating experiments to heavily 
doped n-type silicon samples. We find excellent agree- 
ment between the theoretical and the experimental e-ph 
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FIG. 1: Schematic illustration of the constant energy ellip- 
soids of Si conduction band valleys. The valleys are located 
close to X-point in the first Brillouin zone. Elastic scattering 
rates 1/n and 1/t2 couple the different classes of the valleys. 
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temperature responses. 

As discussed above the electron-phonon coupling in 
semiconductors can be described through deformation 
potential coupling constants, which do not depend on 
the electron variables (in a single valley). The strain in- 
duced conduction band energy shifts Svi (I = 1, 2, . . . , L, 
L is the number of valleys) can be written conveniently 
in matrix notation as Sv = De, where {Sv} l = Svi 
and D is the deformation potential L x 6 matrix (con- 
taining the deformation potential coupling constants), 
e = [ e xx e yy e zz e xy e xz e yz ] is the strain compo- 
nent vector and e a p = i(du a /da + dup/d/3) are the 
symmetric strain components of displacement u. For ex- 
ample, for the six Si conduction band minima (see Fig. 
[TJ wc have Svi = E d (e xx +e yy +e zz ) + E u eu [3, where E d 
(S„) is the dilatational (uniaxial) deformation potential 
constant. 

Here we deal with long wavelength limit where the 
phonon field can be identified with a classical acoustic 
wave u = euexp[— i(q-r — cut)] with polarization e (|e| = 
1). The strain now reduces to t a p = -y-($/3e Q + q a efi)qu 
(q = q/ \q\), which can be expressed in matrix form as 
e = —iqSu and we find equation 



Sv 



-iqDSu, 



(1) 



which describes how the displacement makes the band 
edges oscillate in a many valley system. In the diffu- 
sive long wavelength limit the phonon momentum itself 
cannot transfer the electrons from one minima to an- 
other, because this process would require large momen- 
tum q ~ l/o (a is the lattice constant). Then the electron 
nonequilibrium, generated by the acoustic field, relaxes 
towards local equilibrium by two processes: diffusion and 
elastic intervalley impurity scattering. When the strain 
lifts the valley degeneracy elastic intervalley scattering 
provides a path for the electron system to relax towards 
local equilibrium. This path is favorable if the time scale 
related to diffusion over length ~ g _1 is sufficiently large, 
i.e., when q 2 D < l/r iv , where l/r iv is the total elastic 
intervalley transition rate (D is the diffusion coefficient). 
In this limit the linearized many-valley relaxation-time 
Boltzmann equation ^| reduces to a simple rate equa- 
tion, which couples the change in the electron density Sni 

of valley I to that of valley m via intervalley scattering 
_i 



rate r, 



dSni 



~l [5m - Sn m - 2v 1 (e F )(Svi - Sv m )] 



(2) 

Here V\{ep) is the single spin and valley density of states 
at Fermi level e F ~> ksT. We assume that strain equiv- 
alent valleys are coupled with rate t^ 1 = t^ 1 and that 
the valleys whose degeneracy can be lifted with strain are 
coupled with rate t^ 1 = . In the case of Si the coaxial 
valleys are always equivalent whereas the degeneracy of 
the perpendicular valleys can be lifted (see Fig. . Now 



the solution of Eq. J5J is 5n = 2v±{e F )(l+iuiTi v ) 1 M5v, 
where r iv = t 2 /L and {M} l m = S^ m - L^ 1 . The 
phonon-electron relaxation rate 1 /r q is related to the dis- 
sipated heat Q of the acoustic field through standard re- 
lation l/r q = Q/J £ = u> (lm{Sv} ■ Re{<5n}} /J e , where 
J E is the acoustic energy flux density and ( ) stands for 
time average. Using this relation and Eq. wc find 



(r g xY 



2 Vl {e F ) 



pdv\ 



1 + u 2 qX rl 



(3) 



where we have used linear dispersion relations to q \ — 
v\q (A is the mode index). The factor <&\ — 
e T S T D T MDSe and it obviously depends only on the 
polarization e, on the direction of propagation q and 
on the deformation potential coupling constants. In the 
case of Si we have $a = 22^ [J^ii'l^i) 2 — • e) 2 ] and 
= 6/r 2 . Note that Eq. 10 does not depend on 
screening, because there are no total electron density fluc- 
tuations, i.e., $ n i = 0- 

We can describe a degenerate electron system by an 
equilibrium distribution at temperature T e . This holds 
even in the presence of net heat flow between elec- 
trons and phonons. The heat flow only creates a non- 
equilibrium between the electrons and phonons, which re- 
laxes towards equilibrium at rate 1 / T q \ per single phonon 
mode. By following Perrin and Budd [2(j this non- 
equilibrium can be expressed using the relaxation time 
approximation of the phonon-electron collision integral 



dN(to qX ) 
dt 



ph - 



N{Lu qX )-N Te {u q \) 



(4) 



where N(uj q \) and N T (uj q x) = [1 - exp(hhj q x/k B T)] 1 
are the nonequilibrium and equilibrium phonon distribu- 
tion functions, respectively. The total stationary heat 
flow P through the coupled electron-phonon system is 
the energy average of the collision integral: 



P = 



dq (dN{w qX ) 



(5) 



ph- 



where the summation is performed over the acoustic 
eigenmodes of the crystal. The only experimentally 
meaningful situation is such that the phonon system is 
coupled to some thermalizing bath, which is at temper- 
ature To. If the coupling is strong or P is small we can 
approximate N(u> q \) ~ Nt h (u q \), where T p h is the (pos- 
sibly local) phonon temperature, and Eq. 10 reduces to 
the familiar form: 



P = F(T e ) - F{T ph ), 



(6) 



where F(T) is the energy flow rate control function. 
Using Eqs. I0-© and assuming that (/csT/ftr" 1 ) 2 is 
clearly below unity the energy flow rate control function 
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can be expressed in a closed form 
vi (e F )B 5 



F(T) = 



TT 2 pdfl 5 

2 Vl {e F )~lB b 



{k B Tf 
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r 2 (k B T)\ (7) 



where the first equality is valid for arbitrary many- valley 
system. The constant B 5 — dxx 5 / [1 — exp(x)] = 
1207r 6 /945 and ( ) n stands for average over a solid an- 
gle. The second equality applies for silicon and there we 
have further assumed that the phonon eigenmodes are 
isotropic and that they are described by the longitudinal 
and transversal sound velocities vl and vt- 

Eq. is valid when (k B TjhT^) 2 < 1 and q\D < 
1/t,„. At low temperatures the dominating condition is 
the latter and can be written also as Qt^x/ 1 "™/ 7 " < lj 
where r is the momentum relaxation time. Condition 
qrl^J Ti v /r = 1 defines the crossover temperature be- 
low which elastic intervalley scattering induced electron- 
phonon relaxation dominates over diffusion. If Tj„ is not 
orders of magnitude larger than r this differs very little 
from the impure-pure threshold q?l = 1. 

Eq. suggests that intervalley scattering induced 
electron-phonon energy relaxation rate r" 1 oc Ti V T^ , 
which can be seen from approximate rate equation 
dPj dT e » C e r c _1 , where C e = ^T e is the electron heat ca- 
pacity. As the phonon mediated dephasing rate 1 / rf oc 
rr 1 |ld|| we find an important relation l/if h oc n v T^. 

As already pointed above screening plays no role in 
l/ r qA and as a result intervalley scattering induced 
electron-phonon energy flow rate in Eq. does not in- 
clude any screening parameters, like for example screen- 
ing wave vector k. Note, however, that there exists also 
single- valley contribution to the energy relaxation which 
is due to number density fluctuations in a single valley, 
but this conribution is strongly screened in doped semi- 
conductors [lij. By using the single valley result calcu- 
lated by Sergeev et al. y| and Eq. we find that the 
ratio between many-valley and single valley energy flow 
rate scales roug hly as - 1500 (Ik) (t iv /t) (k/T) 2 , where 
[k] = nm _1 and [T]=K. Thus the many- valley effect is 
expected to fully dominate in the diffusive limit at high 
electron densities and low temperatures. We have tested 
Eq. experimentally in the case of n + Si: 

The n + Si samples were fabricated on unibond silicon- 
on-insulator substrates. Properties of the samples are 
listed in Table|Uand a detailed description about the sam- 
ple fabrication can be found in |2l| . The sample geometry 
and the experiment is depicted in Fig. HJa). In the ex- 
periments the samples were mounted on a sample holder 
of a dilution refrigerator. The electron and phonon tem- 
peratures were simultaneously measured by utilizing the 
superconductor-semiconductor-superconductor (S-Sm-S) 
thermometry |22j while the electron gas in the Si film 



TABLE I: The characteristics of the samples: N - carrier 
concentration, p e - 1.5 K resistivity, I - electron mean free 
path, d - n + Si film thickness. All samples have 400 nm thick 
buried oxide layer. 



sample 


TV (10 25 m~ 3 ) 


p e (10~ 5 fl m) 


I (nm) 


d (nm) 


A 


3.5 


1.04 


5.06 


70 


F 


6.7 


0.63 


5.42 


58 


G 


12.0 


0.51 


4.54 


58 


H 


16.0 


0.44 


4.34 


58 



was heated with a DC power density P — p e J 2 created 
by electric current density J. Note that as the electronic 
coupling to the n + Si film is made via superconducting Al 
the heat flow in the experiment follows accurately a path 
electrons^phonons^substrate/sample holder (phonons) 
and, therefore, the experimental P is equal to the left- 
hand-side of Eq. 0. Heating of the electron gas can 
cause a substantial increase in the temperature of the 
phonon thermometer, as reported recently for a similar 
n+ Si sample as discussed here 0. To assure that the 
nonequilibrium phonon distribution (of the phonons that 
interact with the electrons in the Si layer) can be reason- 
ably described with an equilibrium distribution function 
we consider heating power range where (T p h — Tq)/Tq is 
clearly below unity. 

Fig. Olb) shows the experimental power density vs. 



T a thermometer 
Al 



T„ h thermometer 
J 60 urn 




T 6 -T 6 (10- 3 K 6 ) 

e ph v ' 



FIG. 2: (a) Schematic illustration of the sample geometry 
and the measurement setup. The ~9500 fim long n + Si film 
is heated with a DC current density J. T e and T p h are mea- 
sured using current biased S-Sm-S (Al-Si-Al) contacts (only 
the biasing circuit for T e is depicted). T v h thermometer is 
electrically isolated from the main Si film by a ~1 (im gap. 
(b) The power density P — p e J 2 vs. T® — T® h for samples 
with different carrier concentrations at bath temperature of 
265 mK. 



4 



7.0 


N 


S 


t 2 


1 i 


it , " 


„ 6.0 


■ 3.5 


2.2 


1.9 








6.7 


3.1 


2.2 






r n 
^- o.u 


12 


4.6 


2.7 






"e 


16 


6 7 


3 5 






----- 4.0 




















♦ 




CO 

? 3.0 


















♦ 




















*/ 










2.0 








■ i 






3 4 


5 


6 7 


3 910 


2 



N (10 25 m 3 ) 

FIG. 3: Slopes 5* of the linear fits in Fig. fflb) and interval- 
ley scattering time T2 [determined from S and Eq. J7J] as a 
function of electron density N. The dashed curve is a poly- 
nomial fit that serves as a guide for the eye. The inset shows 
tabulated values of S and T2 (in the units of the axes). 

T e 6 - T£ h at bath temperature T » 265 mK. The solid 
curves are least of square fits to P = S(T® — T® h ) with 
the slope S as a single fitting parameter. We observe that 
the electron-phonon temperature response predicted by 
Eq. J7J describes all the samples extremely well. The 
slopes S are plotted against the electron density in Fig. 
|3I (left vertical axis). S increases as function of N, which 
is expected result, because F(T)/T 6 oc v\(£p). 

In order to perform more quantitative comparison be- 
tween the theory and experiment we estimate the den- 
sity of states from free electron gas expression v\{ef) = 
2,° (ep ) = (m de /2nh 2 ) (3tt 2 N/L) 1/3 , where we use Si 
density of states mass rride — 0.322m e (m e is the free 
electron mass). For the other parameters we use the 
typical values for Si: E u = 9.0 eV, pd — 2.33 x 10 3 
kgm" 3 ,w T(L) = 4700 (9200) m/s. Now the interval- 
ley scattering time T2 = 6r„ can be determined from 
S = F{T)/T 6 [see Eq. 0] and it is plotted on the right 
vertical axis of Fig. [3J The cross-over temperature from 
the condition qxl \/ T{ V / V = 1 is found to be w 5 K (aver- 
age from all the samples). Thus we are at <itI \J Ti V j t <C 1 
limit. 

Eq. 10 gives also the phonon or ultrasonic attenuation 
constant a q \ = (T q \) 1 / (2v\). Using this result and the 
ultrasonic attenuation data obtained by M. Dutoit 0] 
from n + Si with N = 2 x 10 25 m~ 3 at temperature of 
2 K at w = 1.48 x 10 9 1/s we find t 2 (N = 2x 10 25 
m~ 3 ) w 0.3 ps. This fits to our measurements ex- 
tremely well, which is an important result: experiment 
that probes heat transport between electrons and one co- 
herent acoustic mode |l4j coincides with our experiment 
that probes heat transport between electrons and phonon 
gas obeying quantum statistics. 

At high N one would expect slowly decreasing or a 
roughly constant r<i, while our results show a weak in- 



crease as a function of N. This unexpected result could 
be explained by noting that our samples are in the limit 
of strong disorder (kpl w 3.6 on average from Table HJ. 
Whereas, Eq. (JJJ is essentially based on a semiclassical 
free electron gas model, at least finally when the approx- 
imation V\{ef) = v\{ef) is made. Correction terms aris- 
ing from interaction and quantum interference effects can 
be included to our model in the spirit of Ref. As 
similar terms appear in the conductivity the magnitude 
of these quantum corrections can be estimated from low 
field magnetoresistance and temperature dependency of 
resistivity. At the moment such data for thin film n + Si 
is not available. 

Finally, we point out that the intervalley scattering 
induced electron-phonon energy relaxation can be ob- 
served also in several other material systems than n + 
Si. Canonical examples would be n + Ge and two dimen- 
sional electron gas in (111) Si inversion layer. As the 
T-point in the valence band of elemental semiconductors 
is divided into heavy hole, light hole and split-off bands 
the effect should be particularly strong in various hole 
systems. However, due to complicated nature of the va- 
lence band maximum and effectively zero distance of the 
different bands in k-space the theory, which is valid for 
conduction band electrons, should be modified j23l|. 

In summary, we have studied the effect of elastic inter- 
valley transitions on the electron-phonon energy relax- 
ation rate in many- valley semiconductors in the diffusive 
limit. We derived a general expression for the electron- 
phonon energy flow rate [Eq. J7J)] and discussed the spe- 
cial case of n + silicon. Low temperature experiments 
on heavily doped Si samples were performed and good 
agreement between the theory and the experiment was 
found. 
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